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ABSTRACT 

Improved data reduction techniques for 3D data cubes obtained from Fabry-Perot inte- 
gral field spectroscopy are presented. They provide accurate sky emission subtraction 
and adaptive spatial binning and smoothing. They help avoiding the effect analogous 
to the beam smearing, seen in Hi radio data, when strong smoothing is applied to 3D 
data in order to get the most extended signal coverage. The data reduction techniques 
presented in this paper allow one to get the best of both worlds: high spatial resolution 
in high signal-to-noise regions and large spatial coverage in low signal-to-noise regions. 

Key words: methods: data analysis — techniques: 3D spectroscopy. 



1 INTRODUCTION 

In recent years, many Ha kinematic al surveys have 
been done in photon counting mode iHernandez et all 
1200.1 iGach et alJ |2Qq3) using a Fabry-Per ot (FP) inter- 
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may present differences nearly as large as two orders of mag- 
nitude in signal-to-noise ratio (SNR) between starburst re- 
gions and diffuse (inter-arm) regions. Even after hours of 
integration, which means, in fact, about 5 minutes per FP 
channel, the signal in diffuse areas is often too weak to give 
reliable results. 

Since the SNR needed depends on the moments of the 
line profile that has to be analysed (fiux, velocity, dispersion, 
skewness, kurtosis, ...), a typical solution to this problem is 
to apply a spatial weighted averaging (or smoothing) to cor- 
relate neighbouring data points. Doing so, the spatial area 
over which the line profiles can be analysed can be greatly 
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increased at the price of a degraded spatial resolution, hid- 
ing sometimes structures that require high resolution to be 
rendered. Thus, to make it possible to resolve small kine- 
matical details in the galaxy structure and to have the most 
extended coverage possible, it is necessary to create different 
velocity maps with different levels of smoothing. A typical 
approach is to create a non smoothed map, a map smoothed 
with a 3x3 Gaussian kernel and another one convolved with a 
6x6 Gaussian kernel. Then, the data analysis must be done 
three times. Moreover, lowering the spatial resolution in- 
volves a fiux contamination by extending erroneously signal 
from high SNR pixel areas to low SNR pixel zones. 

Given these characteristics. Ha FP observations of 
galaxies would greatly benefit from an adaptive smoothing 
scheme that would be a function of the SNR of the resolu- 
tion element. Simil ar work has been done on integral-field 
spectroscopy (IFS) llCappellari fc Copirl200^ with Voronoi 
tessellation. The main aim of their work was to bin 3D data 
so that the resulting SNR of each bin would tend towards a 
single value. The algorithm they developed has been worked 
on and adapted to scanning FP data, as described in section 
|H] Results are presented in section |1] 

However, while adaptive smoothing can allow increas- 
ing the SNR in weak emission regions without decreasing 
the spatial resolution in high SNR areas, in FP mode other 
aspects of the data reduction must be considered prior to 
smoothing. Section |5| will show that the sky emission spec- 
trum can not be considered homogeneous throughout the 
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Figure 1. Integrated spectra (continuum subtracted) of three 
different sky regions as seen in the data cube of the galaxy 
NGC 5194. Each region is roughly composed of 500 pixels. The 
plots are normalised accordingly. The pixel size on the sky is 
2.56 arcsec'^. The first order of the spectra is between 6569.0 and 
6576.3 A. There are 48 spectral bins for every spectrum. Top ; 
Before sky subtraction. The spectra discrepancies are due to inho- 
mogeneous interference filter transmission. Middle : After sub- 
tracting a median spectrum from the data cube. Bottom : After 
subtracting a sky cube from the data cube. The feature around 
6574. SA comes from the sky emission line at 6596. 64A that is 
not perfectly suppressed. 
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Figure 2. Top : Transmission curves at two different positions on 
the interference filter used to observe NGC 5194 (red and green 
l ines). Th e Ultraviolet and Visual Echelle Spectrograph (UVES; 
IDekker et al. ,,20 00) sky spectrum around the central wavelength 
of the filter is shown as the black line. The spectrum resolu- 
tion have been degraded to the one obtained in the data cube 
of NGC 5194 (R~18600). The FP etalon Free Spectral Ranges 
for this galaxy are also shown. Bottom : The sky spectrum as 
seen in a data cube considering the different transmission curves 
of the filter. The sky spectrum present in every FSR of the FP 
etalon have been folded back into FSRo- 

data cube and an efficient way to subtract it will be pro- 
vided. 



2 SKY EMISSION IN FP DATA CUBES 

2.1 Sky emission inhomogeneities 

Night sky emission lines are very prominent in the red part 
of the visible spectrum. Around Ha at rest, geocoronal Hao 
(6562.8 A) and strong OH sky emission lines (e.g.: 6553.6 
A, 6568.8 A, 6577.2 A, 6577.4 A) contaminate the spec- 
trum. When Ha is redshifted, weak to strong OH lines (e.g.: 
6583.41 A, 6596.64 A, 6603.99 A, 6604.28 A) continue to con- 
taminate the spectrum a t different levels (for a d etailed sky 
emission lines atlas, see lOsterbrock et alJll99^ . Although 
the wavelength and relative intensities of these lines are well 
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documented, it is noticeable that they are not homogenous 
across all the field of view (see top panel of Figure . 

The strength of the OH emission varies throughout an 
observation : 

• The OH radicals' column density changes with the air- 
mass; 

• The OH emission changes as the line of sight crosses 
different OH radicals layers; 

• Radiation (sunrise, sunset, moon) excites the OH radi- 
cals. 

These variations affect large scale OH emission struc- 
tures. Over time, since the FP etalon is scanned at least 
ten times during an observation (and sometimes up to 50 
times), these variations are averaged out. Thus, the main 
contribution to the inhomogeneities of the sky spectrum is 
not the OH emission variation. 

Interference filters are placed in front of the FP etalon 
to restrict the wavelength interval seen by the etalon. These 
filters have a typical profile width of 15 A. Their central 
wavelength is chosen to match the systemic velocity of the 
galaxy. By scanning the filters with a spectrophotometer, it 
has been realised that the filter's transmission shape changes 
spatially. Since these filters are placed in a converging beam 
(rather than into a pupil), their transmission inhomogenities 
are visible in the data cube (Figure |5J. The dissimilarity of 
the sky spectrum observed in the data cubes and the ones 
modelled agree (bottom panel of Figure |21 and top panel of 
Figure 0. Thus, one cannot assume that the sky emission 
spectrum seen in a single FSR of the FP etalon will be ho- 
mogeneous across the field of view as long as the interference 
filters transmission is not perfectly uniform. 

The subtraction of the sky emission is an important 
process. As it will be described in section |5] the kinemati- 
cal information is extracted by looking at only one emission 
line. An improperly removed sky emission will cause inac- 
curate velocities to be found as the amplitude of the sky 
residuals may be greater than the galaxy's signal. In order 
to properly subtract it, one must use a process finer than 
just subtracting a median spectrum (middle panel of Figure 

ID- 

2.2 Sky cube subtraction 

To account for the differences in sky emission seen in the 
data cube, one must construct a sky cube and subtract it 
from the data cube in order to properly remove the sky 
emission. To do so, the data cube must be divided into two 
parts: the galaxy-dominated regions and the sky-dominated 
regions. Then, a sky cube may be constructed by fitting the 
sky-dominated regions and interpolating the sky spectrum 
in the galaxy-dominated regions. The most challenging part 
of the problem is to disentangle the sky-dominated regions 
from the galaxy-dominated ones. 

The galaxy's emission that is redshifted by rotation 
moves spectrally across the data cube. Thus, the spectrum 
that is the most constant across the data cube is the sky 
spectrum. By taking the spectrum of the data cube (the 




Figure 3. Spectrally summed data cube of NGC 5194 observed 
in HcK light. The lighter sections of the image are those considered 
as sky-dominated by the algorithm described in Section l2.2l The 
boxed blue regions are those whose spectra are shown in Fig- 
ure Q Note that the galaxy's companion, NGC 5195, lacks Ha 
emission and has not been considered as galaxy-dominated. The 
bright halo on the left side of NGC 5194 is a ghost coming from 
the bright core of the galaxy. Since this ghost contains emission 
lines from the galaxy, some parts of it are considered as galaxy- 
dominated. Ghost subtraction is a tedious task and it is being 
worked on. 



median of the spectrum that is above the continuum), one 
gets an estimate of the sky emission spectrum structure. 
Then, the data cube is spatially binned to get a fairly good 
SNR on the sky emission. An optimised target SNR was de- 
termined to be ~ 40. Depending on the strength of the night 
sky emission, the total integration time and the pixel size 
on the sky, the size of the bins ranges from ~ 250 to ~ 625 
square pixels, giving an angular size on the sky ranging from 
~ 100 to ~ 1000 arcsec^. Then, the spectrum of these bins is 
cross-correlated with the median spectrum of the data cube. 
However, as discussed in section ITTI the sky spectrum seen 
in the data cube changes from one region to another. Thus, 
in order to determine the sky-dominated regions, one must 
chose a correlation factor that is lower than 1.0. Since the 
apparent wavelength shift of the spectra is typically smaller 
than a channel (~0.2A), a correlation factor of 0.9 gives a 
good approximation of the sky-dominated regions. To bet- 
ter disentangle the galaxy-dominated regions from the sky- 
dominated ones (in cases where the galaxy emission falls 
directly on a sky emission line), the emission amplitude of 
sky-dominated regions must be within ±50 percent of the 
median spectrum. Once the sky-dominated regions are de- 
termined (see Figure O, a 4"^ degree irregular surface fit 
is made to interpolate the sky spectrum into the galaxy- 
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dominated regions. Every sky-dominated region is given the 
same weight. 

Finally, this sky cube is subtracted from the original 
data cube. The result of this process is shown in the bottom 
panel of Figure 



3 ADAPTIVE BINNING OF 3D DATA CUBES 

When reducing a data cube, the main objective should be 
to get the best possible spatial and spectral resolutions. The 
spatial resolution is generally limited by the seeing if the 
physical sampling of the seeing is sufficient. However, it is 
not as straightforward when considering the spectral resolu- 
tion, which depends on the SNR. While it is an easy task to 
reach the theoretical spectral resolution in bright regions, it 
becomes less obvious in low SNR areas. 

Doing a spatial averaging by correlating neighbouring 
pixels is a well known method to enhance the signal cov- 
erage. However, in many cases, this solution is unsatisfac- 
tory since it degrades the spatial resolution even in the 
brightest emission line regions where it is not necessary, 
hiding structures that require high spatial or spectral res- 
olution to be rendered. Moreover, this smoothing induces 
not only a spatial pollution but also a spectral one as it 
transports the flux from the most intense regions towards 
the weaker emission areas, completely drowning the signal 
in th ese areas. For instance, the beam smearing effect in Hi 
data iSwaters. Madore. fc Trewhella 2000 ') produces an un- 
derestimation of the Hi rotation velocity at small galactic 
radii where the velocity gradient is the greatest. 

Hof emission intensities in galaxies display a very large 
dynamic range (several orders of magnitude), from the 
brightest compact and discrete Hn regions to the weaker 
diffuse Ha emission in the inter-arm regions. The higher 
spatial and spectral resolutions can usually be reached from 
the Hii regions without any smoothing. On the other hand, 
the contrast may be very sharp between a bright Hn region 
and a weak diffuse emission area. 

It is thus relevant to introduce the concept of spatial 
smoothing adapted to a given SNR on the same data cube 
allowing maximum spatial resolution in regions of strong 
emission (e.g. Hn regions) and maximum signal coverage 
in regions of weak emission (e.g. inter-arm regions). S uch a n 
approach was recently used by Cappellari and Copin i2003l) . 
working on IFS 2D data. The aim of their work was to bin 
radial velocity maps according the spectra SNR. According 
to them, the ideal binning scheme should give bins that sat- 
isfy the following requirements: 

• Topological requirement: The bins should properly 
cover the sky region, without holes or overlaps. 

• Morphological requirement: The bins should be as 
round as possible to optimise the spatial resolution. 

• Uniformity requirement: The resulting SNR of every bin 
should tend towards a single value. A minimum value should 
be required and no more attempts should be done to accrete 
pixels on a bin that already have the chosen SNR. 

The uniformity requirement is easily fulfilled when the 



resulting SNR of a bin is computed by mean of 

^ Signali 

S N Rhinned j , 

where Signali and Noiset are the signal and noise value of 
the single ith pixel composing a bin. However, as the dif- 
ferent spectral components of every pixel are not taken into 
account, it is impossible to predict the real SNR of a spec- 
trum after its addition to another one. This binning scheme 
also assumes that reaching a greater SNR is simply a mat- 
ter of adding more pixels to the bin. Nevertheless, as spatial 
resolution is limited, accreting pixels that are too distant 
leads to adding pixels that are not correlated together as 
they come from different Line Of Sight Velocity (LOSV) re- 
gions. The binning would be more accurate if the SNR of a 
bin would be recalculated each time a new pixel is added to 
the bin. As the addition of pixels that come from different 
LOSV regions will not increase the SNR of an emission line, 
this will avoid merging structures that are not related. 



3.1 Finding the emission lines 

The determination of the SNR of the emission lines is the 
key of the binning algorithm. Before any attempt is made 
at determining the SNR of an emission line, we first have to 
find where the emission line is located. 

The boundaries of the emission lines are determined 
by a two-pass algorithm. First, the mean signal level of the 
spectrum is found. Then, the highest channel of the spec- 
trum is taken as the starting point of the emission line. If 
the value of the two channels on both sides of the maximum 
are greater than 75 per cent of the mean spectrum value, 
this point is considered to belong to the emission line (the 
threshold value of 75 per cent is chosen to match the am- 
plitude of the two channels each side of the maximum of an 
Airy function that is properly (Nyquist) sampled). If not, 
the algorithm will skip to the next maximum of the spec- 
trum until one is found that fulfil the restriction. This max- 
imum is the starting point. This mechanism has been put in 
place when processing CCD data. A cosmic ray could have 
been the maximum point and it must be avoided to look for 
an emission line around it. In photon counting, where the 
observation is multiplexed and the camera is operated at a 
high frame rate, this can not happen. This mechanism also 
helps finding emission lines in really low SNR spectra. 

Then, by walking both sides of the maximum, the first 
points that pass below the median flux and have a change in 
the sign of the slope are found. This delimitates the rough 
emission line. If a change in the sign of the slope is never 
encountered, the first point that falls below the median is 
taken as the boundary. Using those boundaries, the contin- 
uum is computed by taking the mean value of the channels 
that fall outside the emission line. The dispersion of that 
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Figure 4. Comparison of the emission line position determination algorithm (plain curve) and a Gaussian fitting algorithm (dotted 
curve). The signal dispersion relative to the emission line height is plotted on the x axis. The lower line is the mean absolute position 
error (in FP channel) and the upper line is the mean absolute position error plus the position error dispersion (Icr). Left : Results 
without sky residuals. Right : Results with sky residual. Sky residuals were 0.2 the amplitude of the emission line. 



first continuum is computed as follows 



n - 1 



The second pass involves starting back from the maxi- 
mum channel previously found and walking both sides of it 
to find the first points that pass below the continuum plus 
dispersion, but if a change in the sign of the slope is encoun- 
tered in the 2g range of the continuum, this point is taken 
as the boundary. This allows a better avoidance of the sky 
residuals that could still pollute the galaxy's spectra. The 
boundaries found are taken as the true boundaries of the 
emission line. 

Now that the emission line is isolated, it is possible 
to extract the parameters of the spectrum. The continuum 
is computed by taking the mean signal outside the emis- 
sion line and its dispersion is calculated. The centre of the 
emission line is simply computed as the intensity weighted 
centroid of the spectral bins falling into the emission line 
boundaries, with the continuum subtracted. The amount of 
spectral bins used depends on the width of the emission line. 

Figure 0] shows how this selective intensity weighted 
mean algorithm performs as compared to a Gaussian fitting 
algorithm (gaussfit in IDL). The simulated spectra were 
made of a Gaussian convolved with an Airy function. Noise 
with a normal distribution was then added to the spectra 
to account for the many kinds of noise that can be encoun- 
tered. For spectra where the signal dispersion relative to the 
amplitude of the emission line was > 0.2, the Gaussian fit- 
ting algorithm diverges completely. For less noisy spectra 
(0.004 < Oreiative < 0.2), the Gaussiau fitting algorithm 
performs slightly better. However, the additional error in- 
duced by the selective intensity weighted mean is less than 
0.15 channel. 

In the right panel of Figure |1] the same comparison 
is made but it uses a spectrum that is polluted with a sky 



residual that has an amplitude that is 20 per cent that of the 
emission line. As expected, the Gaussian fitting algorithm 
performs better for low noise spectra (0.01 < a-j-eiative < 
0.1), with an absolute error that is 0.2 channel lower. How- 
ever, the Gaussian fitting algorithm is still confused by noisy 
spectra and diverges. 

It was decided to use the selective intensity weighted 
mean algorithm since it performs very well for noisy spectra. 
This signal regime is usually seen in very diffuse areas and 
the ability to resolve kinematics in these regions dictates 
how extended is the signal coverage. 



3.2 Noise sources and SNR 

Due to the absence of read-out noise, FP observations at Ha 
in photon counting mode are dominated by photon noise. 
One component of the SNR must thus be the noise of the 
sum of the photons (N) contained in the emission line that 
are above the continuum (the photons that create the emis- 
sion line itself). Also, to account for the improperly removed 
sky lines, it is assumed that the continuum should be flat 
and its dispersion (a) is used as another source of noise. This 
dispersion can also account for the read-out noise of the de- 
tector while processing CCD data. Obviously, in this case, 
the CCD signal must be gain-corrected so that the data is 
expressed in photons. The SNR is computed as follows 



SNR ■ 



N 



where 



is the number of photons in the emission line and c is the 
continuum level. 



© 2006 RAS, MNRAS OOO.ITIITol 



6 0. Daigle et al. 




-3M5^56^ 15^54^ 15^52^ 15^50= 15^48= -5-^46^ '5^44^ 13M5"^56= 15^54= 15^52^ 15^50^ 15'"4Ss -5-^46^ '5^44= 

RA (J2000} RA (J2000) 



Figure 5. Comparison of the 2D rotation velocity maps of NGC 5055 at the centre of the galaxy with different smoothing algorithms. 
Top left : No smoothing. Top right : Gaussian 6x6 kernel. Bottom left : Adaptive smoothing, target SNR of 7. Bottom right : Bin 
sizes obtained with the adaptive smoothing. 
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Figure 6. Comparison of the full 2D rotation velocity maps of NGC 5055 with different smoothing algorithms. Top left : No smoothing. 
Top right : Gaussian 6x6 kernel. Bottom left : Adaptive smoothing, target SNR of 7. Bottom right : Bin sizes obtained with the 
adaptive smoothing. 
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3.3 Binning and smoothing 

The binni ng alRorithm used is an imp rovement of the one 
defined bv ICaopeUari fc Copinl 1)200,^ . with the difference 
here that it allows the SNR of the bin to be recalculated 
every time a new spectrum is added. This method leads to 
a greater scatter in the resulting SNR of the bins, but, as 
previously stated, it allows the limitation of the spatial reso- 
lution to be taken into account. This means that the SNR of 
the bins will not necessarily grow as more spectra are added. 
It was found that targeted SNR values in the range of 5-10 
are optimal (SNR values higher than 10 creating bins that 
are too large and into which the SNR is not increasing as 
more pixels are added). The original algorithm was also op- 
timised to allow the processing of large data cubes (typically 
512x512 pixels and 48 channels) in a reasonable amount of 
time (e.g. less than 10 minutes on a typical personal com- 
puter) . 

After the bins are delimited, to avoid singularities at the 
bins boundaries, a Delaunay triangulation algorithm is used 
to smooth the edged of the bins, taking the fiux weighted 
centroid of the bins as the nodal points of the triangles. 
This smoothed data cube is finally used to extract the radial 
velocity map, using the algorithm described in section 13.11 
to find the barycentre of the emission line. 

Applying an accurate sky emission subtraction is very 
important before attempting to adaptively smooth a data 
cube. As the SNR of a spectrum is calculated by looking at 
only one emission line, an inadequately subtracted sky will 
cause the algorithm to trigger on the sky residual rather 
than on the galaxy's signal. The resulting smoothed cube 
will therefore have a constant SNR on the sky residual and 
the SNR of the galaxy's emission will rarely be optimum. 



4 RESULTS 

Figures|S]and|S]show the results of different levels of smooth- 
ing on the extracted velocity maps of the galaxy NGC 5055, 
using Ha. Figure 13 concentrates on the central parts of 
the galaxy, where the accuracy of the kinematical informa- 
tion is crucial (see section 14. H . The high resolution (non 
smoothed) velocity map gives accurate kinematical informa- 
tion but has a very small spatial coverage. When a Gaus- 
sian smoothing is applied, these small kinematical details 
are lost, rendering the radial velocity map inaccurate. Note 
that the unsmoothed map has been thresholded using its as- 
sociated monochromatic map to avoid having pixels of low 
monochromatic flux, whose velocities were incorrectly de- 
termined, polluting the map. A threshold of 5 photons has 
been used. 

When considering the extended part of the galaxy (Fig- 
ure ISJ, it is obvious that spatial smoothing must be ap- 
plied in order to resolve the kinematical information at large 
galactocentric radii. In that strong smoothing per- 

mits to have the best signal coverage but still causes velocity 
inaccuracies in the high SNR regions. 

In both cases, the adaptively smoothed map gives the 
best results: high resolution, high level of detail, great cov- 



erage at the centre of the galaxy and wide signal coverage 
in the low SNR regions. The median residuals of the adap- 
tively smoothed map (Figure are lower than the ones ob- 
tained with other smoothing algorithms (7.7 kms~^ against 
8.9 kms~^ and 10.7 kms^^for the Gaussian 3x3 and 6x6 
smoothing) . 

Naturally, this method is efficient when the SNR varies 
greatly from one region to another. If the SNR is the same 
everywhere, the effect of the adaptive smoothing will be 
roughly the same as a fixed size kernel smoothing. 



4.1 Effect on rotation curves 

Rotation curves (RC) derived from the data cubes are shown 
in Figurcs|Hland|5| The RC were extracted using the rotcur 
(Bcecman 1987) routine of the GIPSY package, using fixed 
inclination (i) and position angles (PA) at all radii. The er- 
ror bars are the rotation velocity errors as determined by 
rotcur. When there is little kinematical information into a 
ring or when there is a great scatter in the velocities com- 
posing the ring, the errors are large. The deep notch in the 
non smoothed rotation curve around 10" is explained by the 
few emission points in the ring at that radius, as shown in 
the bottom panel of the figure, and is thus inaccurate. When 
a slight smoothing is applied, this feature disappears. How- 
ever, when the smoothing is too strong as it is the case for a 
6x6 Gaussian kernel, the rising part of the rotation curve is 
lowered. This produces an effect similar to the beam smear- 
ing observed in Hi data due to the limited spatial resolu- 
tion. Adaptive smoothing should thus be preferred since, as 
shown bv .Blais- OucUcttc ct al. ( 1999), the rising part of the 
rotation curve is crucial in the determination of the dark- 
to-luminous mass ratio and of the mass models parameters. 
Note that the rotation curve extracted from the adaptively 
smoothed radial velocity map is very close to the rotation 
curve extracted from the unsmoothed one within their error 
range (±1(t) for nearly all points. 

When considering the extended parts of the rotation 
curve, strong smoothing must be applied to resolve the kine- 
matics, as the signal in these regions is very weak. Figure 
|U] shows that very little kinematical information is avail- 
able when no smoothing is applied. Thus, when considering 
fixed-size kernel smoothing, multiple rotation curves must 
be merged together to recover both resolution and signal 
coverage. Also, the Gaussian smoothing sometimes makes 
small detail of the RC disappear by convolving uncorrelated 
data points (like the structure around the galactocentric ra- 
dius of 150" which corresponds to a spiral arm of NGC 5055). 
The fall off seen for the adaptively smoothed curve at radii 
> 250"can be explained if the disk of NGC 5055 is warped 
jBosmalligSlh . This has not been modelled here. 

The advantage of using an adaptive binning/smoothing 
algorithm is obvious for rotation curves extraction. It is the 
algorithm that produces the most accurate rotation curves, 
based on the largest spatial coverage possible. Moreover, it 
allows rendering high resolution kinematical information at 
the very centre of the galaxies where even an unsmoothed 
radial velocity map could not produce accurate results. Hav- 
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Figure 7. Histograms of the absolute velocity residuals of the ve- 
locity maps (smoothed — non-smoothed) of NGC 5055 obtained 
with different smoothing algorithms. The vertical lines show 
the median residual for each algorithm. For the Gaussian 6x6 
and 3x3 algorithms, the median residuals are 10.7 kms~^ and 
8.9 kms~^ respectively. The median residual for the adaptive 
smoothing is 7.7 kms~^ . Note how the histogram of the residual 
of the adaptive smoothing peaks towards kms~^ . 



ing a data cube that is smoothed to a constant SNR would 
also allow one to search for multiple emission line compo- 
nents with greater ease. 



5 CONCLUSION 

It has been shown that rigourous data reduction techniques 
applied to 3D FP data help getting as much information as 
possible out of data cubes with the best possible accuracy. 
Whilst good sky emission removal methods are necessary to 
avoid contaminating the weak galaxy's signal, a smoothing 
algorithm that adapts itself to the signal available helps pre- 
serving spatial resolution and optimises signal coverage. The 
effects of these methods are of high scientific importance: 
they help providing accurate and extended kinematical in- 
formation, as demonstrated in section I^TTI 

It is clear when comparing the three smoothing schemes 
in Figures |S| and |S| that the adaptive smoothing is doing as 
well as the heavy (6x6) Gaussian smoothing in the regions 
of weak signals while preserving the kinematical details in 
the regions of strong signals. Moreover, this technique avoids 
underestimating the velocities (beam smearing effect) in the 
inner parts of the rotation curves where the velocity gradient 
is strong, thus providing much more accurate kinematical 
data. 

These data reduction techniques have already been 
applied for data presented in survey observation papers, 
such as the Big Hq kinematics sample of BARred spiral 
galaxies (Hernandez ct a,l. 200,5) , the Hq Kinemati cs of the 
SINGS Nearby Galaxies Survey llDaigle et alJ 2006h and the 
Virgo high-resolution Ha kinematical survey i Chemin et alJ 
|200y). The data for these surveys were gathered at the Ob- 
servatoire du mont Megantic 1.6-m telescope, the Canada- 
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Figure 8. Comparison of different smoothing algorithms applied 
to a galaxy velocity map and the result on the inner part of the 
extracted rotation curve. Top : The rotation curves. Note that 
the adaptive smoothing algorithm minimizes the beam smearing 
effect. Bottom : The number of points of the velocity map used 
in each crown of the tilted ring model to extract the rotation 
curve. Note that even at the centre of the galaxy, very few points 
are available in the non smoothed map. 



France-Hawaii 3.6-m telescope and the ESO/La Silla 3.6-m 
telescope. Not only do these methods provide data of high 
quality, they also reduce the time spent processing the data 
itself. In fact, a great deal of time was needed to manu- 
ally clean up radial velocity maps where the sky emission 
had been improperly removed. With the growing number of 
galaxies observed in Hq surveys, these techniques helps the 
processing of extensive amounts of data in a short period of 
time. 

The data reduction techniques presented in this paper 
focus mostly on FP data obtained with an Image Photon 
Counting System at Hq. However, these methods can also 
be applied to other data sets, in particular to any 3D integral 
field data obtained with spectrometers based on grisms sys- 
tems (optical fibres, microlenses, slicers) . All the techniques 
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Figure 9. Comparison of different smoothing algorithms appUed to a galaxy velocity map and the result on the whole extent of the 
extracted rotation curve. Points beyond 200" were removed for the non smoothed rotation curve as they were not accurate and were 
distracting. Top : The rotation curves. Bottom : The number of points of the velocity map used in each crown of the tilted ring model 
to extract the rotation curve. 

is provided so that the compilation of the module un- 
der different flavours of Unix operating systems should 
be made easy (including Linux and Mac OS X). The 
routines are also working in IDL under Windows. The 
routines can use as input files raw observational data 
(usually ADA files, containing only photon positions for 
the observation and ADS files for the calibrations. See 
http://www-obs.cnrs-mrs.fr/adhoc/adhoc.html for AD* 
files format.) or they can be feed with interferograms cubes 
that have not yet been phase-corrected, phase-corrected 
cubes, phase-corrected and sky-subtracted cubes, and so on. 
This makes it possible to process data using another data 
reduction package and use these routines for some phases of 
the data reduction only. Other file types (such as the more 
standard FITS format) can be supported with the supplied 
FITS to AD* and AD* to FITS converters. The software 
also provides mosaicing and merging facilities to allow for 
large galaxies and galaxies observed at different times to 



described could also be applied to CCD FP data cubes. 
The sky emission subtraction method described in section 
|5| should work for any data cube in which sky-dominated 
regions are apparent. In the case where the galaxy occupies 
the whole field of view, an observation of the sky must be 
interlaced with the observation of the galaxy and the sky 
cube may then be subtracted from the galaxy cube. 



APPENDIX I: AVAILABILITY 

The data reduction routines described 
throughout this paper are available at 
^http : //www. astro .umontreal . ca/f antomm/reduction 
in IDL language. The intensity weighted mean algorithm 
has been written in C language for performance reasons. 
The C module is loaded in IDL through the Dynami- 
cally Lodable Modules (DLM) mechanism. A Makefile 
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Table 1. Observational parameters 



Parameter 


NGC 5055 


NGC 5194 




q A /rslbr 
Ort\_l o ) uc 


O -£1. 1 o t U lJt;t, 


V sys V / 


504 


463 




14/0S/2004 


J-O/ VJ'J / iiVJUtJ 






1 3h30m07 3s 


Field DEC ♦ i f ,120001 


+42d00mll.3s 


+47dl2mll.0s 


Field of view (') 


13.8 X 13.8 


13.8 X 13.8 


Pixel size (") 


1.61 X 1.61 


1.61 X 1.61 


Focal ratio (F/D) 


2.0 


2.0 


Filter \central (A) 


6584 


6581 


Filter FWHM (A) 


15.1 


19.8 


Filter transmission (%) 


74 


60 


FP interference order (p) 


765 


899 


FP FSR (kms-i ) 


391.88 


333.47 


Number of channels 


48 


48 


Channel step (A) 


0.15 


0.18 


Calibration Finesse 


16.8 


20.7 


Resolution (^) 


12852 


18609 


Total exposure (min) 


128 


246 


Exposure per channel (min) 


2.66 


5.14 



be processed. This package does not (yet) include visual- 
isation routines or any kind of interactive user interface. 
Visualisation applications of the ADHOCw package can be 
used (http : //www-obs . cnrs-mrs . f r/adhoc/adhoc .html ) 
to view AD* files. They are designed to run under the 
Windows operating system, but they have been reported 
to work fine under Linux using the Windows Emulator 
(WINE) package jhttp : / /www . winehq . comTt . 

APPENDIX II: OBSERVATIONAL 
PARAMETERS 

The observati onal data used in this paper were gath ered 
with FaNTOmM jHernandez et al.ll200 A iGach et al.ll2002^ . an 

integral field photon counting spectrometer using a Fabry- 
Perot etalon. The observational parameters for galaxy used 
as case studies in this paper are presented in Table Q 
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